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Abstract 

We report a theoretical analysis of the phonon thermal conductance, k(T), 
for single wall carbon nanotubes (SWCN). In a range of low temperatues 
up to WOK, k(T) of perfect SWCN is found to increase with temperature, 
approximately, in a parabolic fashion. This is qualitatively consistent with 
recent experimental measurements where the tube-tube interactions are neg- 
ligibly weak. When the carbon-carbon bond length is slightly varied, k(T) is 
found to be qualitatively unaltered which implies that the anharmonic effect 
does not change the qualitative behavior of n(T). 
PACS numbers: 61.46.+w,44.10.+i,63.22.+m 



1 



I. INTRODUCTION 



Since its original discovery, carbon nanotubes [1] have received a great deal of attention 
due to fundamental physical interest on nano-scale systems, as well as due to nanotubes' 
potential for useful industrial applications [2]. A very important recent advance has been 
the fabrication of high-purity crystalline bundles of nearly mono-disperse single wall carbon 
nanotubes [3-5]. This allows better experimental control and produces accurate data. It 
also provides opportunities and points to new directions for theoretical analysis of SWNT. 
So far, the electronic and mechanical properties of carbon nanotubes have been extensively 
investigated while there also exist several new measurements of thermal properties of these 
systems [6-9] . The purpose of this work is to present our theoretical and numerical analysis 
of thermal conductance of SWNT. 

Of particular interest to this work is the recent experimental measurements of thermal 
conductance of nanotubes k(T) [6]. This quantity describes the thermal current induced 
by a temperature gradient, and we will analyze a two-terminal measurement of it. Hence, 
maintaining temperatures at the left and right lead to be T L , T R respectively, k(T) is defined 
as k(T) = Q/VT, where Q is the thermal current flowing through the nanotube. The 
experimentally measured thermal conductance of carbon nanotubes indicates that the most 
essential contribution comes from phonons [6,7]. The data of Ref. [6] on aligned multiwall 
nanotubes suggested that at low temperature up to ~ 120K, n(T) can be well fit by quadratic 
form in T, i.e. k(T) ~ T 2 . On the other hand, the measurements [8] of specific heat C(T), 
which is proportional to k(T), also showed a power law dependence, however with a power 
larger than unity but slightly less than two [8]. Finally, measurements on crystalline ropes 
of SWNT indicates [7] a linear temperature dependence up to 30K and an upward bend 
near 30K on the k(T) versus T curve [7]. How to interpret these experimental observations 
presents a challenge to theory. Indeed, although there exist a vast literature on theoretical 
analysis of thermal transport of both electrons and phonons in the context of mesoscopic 
physics [11-19], direct investigation on carbon nanotubes is rare [8,10,20]. 
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The precise temperature power law in thermal conductance of nanotubes is expected to 
depend on, among other things, the detailed phonon dispersion. In this work we will calculate 
the dispersion for several SWNT and investigate its consequence. We neglect electrons as 
their contribution to k(T) of nanotubes can only be observed at temperature less than IK 
[10]. Our analysis of SWNT is based on the Tersoff- Brenner /s empirical potential [21,22] 
for carbon to calculate the phonon dispersion. We then apply a Landauer-Buttiker-type 
formalism to compute the lattice thermal conductance in various SWNTs, for which the 
thermal transmission coefficient is assumed to take a Breit-Wigner form. The predicted 
k(T) shows a quadratic form in temperature for both zigzag and armchair SWNT. Using 
one experimental data to fix an overall shift of k(T), our predicted k(T) is consistent with 
the experimental data of Ref. [6] as the tube-tube interactions are so weak that we could 
compare our calculations for single wall carbon tubes with the experimental results on 
multiwall tubes. As pointed out in Ref. [7], since the tube-tube interactions are also weak, 
the measured k(T) is linear below 30K and shows an upward bend slightly near 30K, i.e. 
the curve appears to be also parabolic-like from 8K to 100K. When the carbon-carbon 
bond length is slightly varied, our results indicate that k(T) does not change qualitatively, 
implying that the anharmonic effect does not alter its qualitative behavior. 

The rest of this paper is organized as follows. In the next section we outline the analysis of 
phonon dispersion and n(T). Section III presents the numerical data while a brief summary 
is given in the last section. 

II. LATTICE THERMAL CONDUCTANCE 

We start by writing down a multi-probe formula for thermal current transmission in 
the same spirit as the familiar Landauer-Buttiker formula [23] for electron transport. As 
discussed above we neglect contributions from electrons and deal with the lattice vibration 
as a phonon gas. The thermal current is therefore driven by a temperature gradient and 
carried by phonons from the probe to another. The thermal current can be written as 
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Q a = \fdEE J2ME)A a p(E) (1) 
where a, (3 label the leads, /i is the Planck constant, E is the phonon energy, and fp{E) = 



\j(e h BT(i _ j s ^he phonon distribution function with temperature Tp at lead (3 . The 
thermal transmission function A a p(E) is defined as 

A afi {E)=S a p-^ aP {E)S a p{E), (2) 

where S a p is the scattering matrix for phonons. Eqs.(l,2) are similar to the Landauer- 
Buttiker formula [23] for electron transport. Let the energy of the s-th branch of phonon be 
hoj s . The thermal current can be rewritten as 

Q a = ^2 / duJs huJs Hf^ huJ s) A a/3 (huj s ) (3) 



or 



I r°° 

Q a = ^T,J d v ■ v *(q) ^-(q) E M<D A «M), (4) 

s f3 

where v s (q) = V q w s (q) is the group velocity. 

For the special case of two probe measurement in one dimension, one can easily confirm 
that Eq.(4) is reduced to the form which has been used to study the quantized thermal 
conductance of dielectric quantum wires [24] . We have 

Q = E/ ^huj s (k)v s (k)(f R -f L )Uk), (5) 

where u> s (k) and v s (k) = 9u1 q^ are the frequency and the velocity of normal mode s of the 
wire with wave vector k, respectively, ( s (k) is the phonon transmission probability through 
the wire, and f a is the phonon distribution function at the right or left lead indicated by 
a (= R,L). We refer interested readers to Refs. [25,26] for more details about the present 
formulation. For the SWNT the thermal conductance, k = Q/VT, can be numerically 
computed through the following form 



K = — ]T / dcuhcu(f R (cu) - f L (cu))C s (u;)/AT. (6) 



where AT = T R —T L is the temperature difference between the two leads, AL is the length of 
the SWNT, and 

Wmax is the maximum of phonon modes in SWCNs, which was investigated 

elsewhere [27]. 

To proceed, it is clear that the phonon frequencies of SWNT at T point and the trans- 
mission function ( s (u) are required. There are a number of methods which can be applied 
to obtain the phonon frequencies, and we will present our calculation of this quantity in 
the next section. On the other hand, it is rather complicated to accurately determine the 
phonon transmission function. For a perfect infinite harmonic lattice the thermal conduc- 
tivity diverges because every mode transmits perfectly. The contacts play an important role 
as in case of electron transport. Even for perfect phonon transmission, the thermal conduc- 
tance is finite due to the contact resistance. Defects, boundary scatterings, phonon-phonon 
interactions and other scattering mechanisms in the conductor inevitably introduce thermal 
resistances. The transmission function ( s (ui) is contributed by lattice imperfection and con- 
ductor boundary, and the anharmonic effect due to interactions can be included in terms 
of the renormalized temperature-dependent dispersion relation. To obtain ( s (uj) for SWNT 
one should in principle solve the scattering of a phonon wave by the atoms of the nanotube 
which is a nontrivial problem not solved so far. However if the length of a conductor is 
much larger than its cross-section size, Cs(w) can be approximately obtained by solving lin- 
ear elastic equations assuming the conductor to be an uniform elastic medium [24]. In ID, 
the continuum mechanics results indicate [24] that ( s (uj) is dominated by sharp resonances 
reminiscent of a resonance transmission of phonon modes: ( s (uj) = 1 at resonant energies 
and sharply becomes much smaller for other energies. Guided by this observation and for 
simplicity of analysis, we shall investigate the consequence of such a resonance transmission 
of phonon modes for carbon nanotubes by parameterizing ( s into a Breit-Wigner form, 



D 2 

CiE)= (E-Eo)* + D*> (7) 
where D and E are the width of the resonance and position of a resonance. Both of these 
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parameters are closely related to the effects of defects and boundary scatterings. Further- 
more, instead of directly considering anharmonic effects, we will investigate its consequence 
by simply varying the carbon-carbon bond length to observe how much the thermal conduc- 
tance can vary. 

III. LATTICE VIBRATIONS 

To calculate the phonon spectra u s at the Y point for armchair and zigzag carbon nan- 
otubes, we shall make use of the Tersoff- Brenner empirical potentia [21,22] for carbon. 
Although this empirical potential only gives an approximation to the vibration frequency 
of nanotubes, it does capture the basic and qualitative characteristics of their phonon spec- 
tra. Tersoff-Brenner's potential is also numerically efficient to allow investigations of large 
systems which is required for our study. The parameters for this potential have been deter- 
mined before and there are two sets of them which have been widely used [22] for structure 
analysis of carbon: 

R$k D$eV Pij/k Sij 5, R^ ] k R$k a c d 

1 st 1.315 6.325 1.5 1.29 0.80469 1.7 2.0 0.011304 19.0 2.5 
2 nd 1.39 6.0 2.1 1.22 0.5 1.7 2.0 2.0813 x 10~ 4 330 3.5 

The meaning of the symbols refers to Ref . [22] . We have checked that for S WNT the second 
set of Brenner/s parameters gave results closer to those obtained by using Tersoff parameter 
[21]. We shall use in our calculations the second set values shown in the table. In addition, 
experimental measurements have shown that phonon modes of SWNT have frequencies up 
to 1700cm -1 . We have numerically confirmed that if the equilibrium distance parameter 
R\f = 1.33 A, which is a value in between the first and second Brenners/s parameter sets 
(see table), the obtained vibration frequencies at T point is consistent with experimental 
data [28] . On the other hand, all the parameter sets give similar total binding energy. For 
instance, for the armchair (4,4) nanotube, the binding energy per atom is found to be -6.8920, 
-7.0881, -6.9993 and -7.0456 eV using Tersoff' s potential, the first Brenner's potential, the 
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second Brenner's potential and our modified second Brenner's potential, respectively. Fig. 
1 presents a comparison between these empirical potentials. In all subsequent numerical 
analysis we shall utilize the revised value of R^f . 

The lattice vibration frequencies of the SWCN at T point can be obtained by diagonal- 
izing the dynamical matrix given by 

D a p{k, Ik') = ^E^/3 (0*. lk ') e^ R W- R ^'\ (10) 

where q is the wave vector, the perturbation of R a (lk) is taken to be O.OOlA, and 

d 2 E 

^ {lkJ ' k,) = dR a {lk)dR p {l'k'Y 

Once the phonon modes are obtained, we are able to calculate thermal conductance by Eq. 
(6). Here we would like to point out that although there have been a number of studies 
dedicated to the lattice vibrations of carbon nanotubes so far [2], it is not suitable for us 
to adopt directly the phonon data in literature. This is because in our situation we need 
frequency data for different sizes of carbon nanotubes as well as the data of frequencies 
when the carbon-carbon bond length is simply varied in order to retrieve the anharmonic 
effects. In principle, we can use ab initio molecular dynamics method to get the frequencies 
if the system is not too big. Since in the present case the system contains hundreds of 
atoms, it is not easy to obtain useful data on the basis of ab initio simulations, and for a 
practical purpose it is better to invoke an empirical potential to compute the phonon spectra 
of SWNCs with larger sizes. 

IV. RESULTS 

As discussed in the introduction, we are investigating long SWCNs such that the left 
end is at temperature Tl and the right end is at temperature Tr. Let us first consider 
the basic characteristic of k(T). At very low temperatures where only the massless modes 
(i.e. u s (0) = 0) are relevant, k/T is simply a constant independent of temperature. At 
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slightly higher temperature, the high frequency modes set in because a SWNT is a quasi- 
one dimensional system in which each unit cell contains so many atoms that there are 
many modes with nonzero cutoff frequencies. Therefore k/T will deviate from the constant 
behavior when temperature is increased. 

The inset of Fig. 2 shows the temperature dependence of thermal conductance for the 
armchair (11,11) and the zigzag (11,0) carbon nanotubes when the transmission is perfect, 
i.e., by setting q s = 1 for all the phonon modes. For these cases our results indicate that k(T) 
increases, approximately, in a parabolic fashion as a function of temperature. We found that 
anharmonic effect has only a slight influence on the temperature dependence of k(T), and 
this is shown by calculating k(T) for different values of carbon-carbon bond length which 
generates different phonon spectra at the T point. Fig. 2 shows k(T) for three different 
bond lengths 1.42A, 1.44A and 1.46A for a zigzag (11,0) system. This result clearly shows 
that the basic shape of k(T) is not changed qualitatively by varying the bond length. The 
magnitude of k(T) is however dependent on bond length as expected: the longer the bond 
length the lower the thermal conductance value. 

In a real SWNT system the flow of thermal current is not perfect due to scattering of 
phonons by various lattice imperfections. It is thus inevitable that <, s < 1 in general and only 
at resonances [24] do we have <, s = 1 [see Eq. (7)]. Nevertheless, as long as the transmission 
resonances are sharp for a large number of modes, which is the case shown by previous 
numerical calculations for various one-dimensional wires [24], we do not expect qualitative 
changes of behavior from that shown in Fig. 2. On the other hand, it is interesting to 
examine the opposite situation, namely the transmission is dominated by a single resonance, 
i.e. the transmission coefficient q s is centered at a single mode E with a resonance width 
D [see Eq.(7)]. Therefore for a small D only very few terms in the summation of Eq.(6) 
have substantial contribution to k{T), while a large D should give a similar result as those 
of Fig. 2. Fig. 3 shows this behavior for a number of values of D by fixing E = D for 
simplicity. Indeed, a large D which essentially gives q s = 1 for a wide range of modes, gives a 
parabolic temperature dependence for k(T). For small values of D , on the other hand, not 
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only the value of k is reduced because only a few modes contribute, but also the temperature 
dependence changes qualitatively. This exercise indicates the importance of the quality of 
phonon wave transmission. 

Recently the thermal conductance of the carbon nanotubes has been experimentally 
measured [6,7]. In both experiments, the results show clearly that the tube-tube interactions 
are quite weak, implying that it would be reasonable at some extent that is also capable of 
comparing the measured results for multi-wall carbon nanotubes with the calculated ones 
for single wall tubes without generating qualitative deviations, as it indeed is. Our result is 
found to be qualitatively consistent with the experimental data [6] as shown in Fig. 4 where 
the circles are the measured data and the dotted and dashed lines are our calculated results. 
The solid line is a least square fit to the experimental data in the range of 10-K" to 120K which 
gives k = 0.00089 * T 1 - 97 (W K^ 1 M~ r ) . Hence the experimental data indicate a parabolic 
dependence on temperature. Surprisingly, our calculated value of k is even quantitatively 
reasonable considering the various approximation used in our analysis: since the Lorentz 
number =jg- - 10~ 12 ( JK^S- 1 ), the size of nanotube is on the order of nm ~ 10 9 (M), 
and at T = lOOfT our calculated k ~ 100 in units of Lorenz number (e.g. Fig. 2), gives the 
same order of magnitudes in magnitude as the measured data. To plot our calculated curves 
in Fig. 4, we have used one experimental data point at T = 120/T to fix the scale of k. It is 
clear that for both zigzag and armchair tubes, the predicted curves qualitatively agree with 
the measured data quite reasonably. In Ref. [7] , the measured k(T) is linear below 30K and 
an upward bend slightly near 30K, i.e. the curve appears also to be parabolic-like from 8K to 
100K. The linear dependence below 30K in Ref. [7] may be the results of particular phonon 
scattering mechanisms or the imperfection phonon transmission. As a result, our results 
appear to suggest that the thermal transport in nanotubes might primarily be mediated by 
lattice vibrations, and the electron contribution is less important. 
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V. SUMMARY 



We have studied the properties of the thermal conductance of single-walled carbon 
nanotubes by means of a Landauer-Buttiker-like thermal conductance formalism using a 
modified Tersoff-Brenner potential. We found that the thermal conductance of the SWCN 
increases almost parabolically with increasing temperature at perfect transmission at low 
temperature. Our results are qualitatively consistent and quantitatively on the same order 
of magnitude as the measured data [6,7] because the tube-tube interactions in these carbon 
nanotubes are weak. This allows us to suggest that in the experimental device the phonon 
transmission is quite ideal, which is consistent with the fact that the high quality of the 
nanotubes [6,7] was ensured in the measurements. It also suggests that it is the phonon con- 
tribution which gave, to a large extent, the observed behavior of temperature dependence 
of the thermal transport in SWCNs. In addition, our results suggest that the anharmonic 
contribution to the thermal conductance of SWCN is only quantitative but not qualitative. 

Although the simple analysis presented in this paper has given a qualitative understand- 
ing of the experimental data of Ref. [6,7], a number of further improvements to theory are 
desirable. An extremely difficult subject is the calculation of phonon transmission func- 
tion for carbon nanotubes. Our choice in this work is the phenomenological Breit-Wigner 
formula motivated by one-dimensional theory of linear elasticity. For realistic nanotube de- 
vices, especially those with lattice imperfections, it will be useful to accurately compute the 
transmission coefficient. Another improvement to the present work involves a more accu- 
rate analysis of the phonon frequency, perhaps from first principle methods. Finally, since 
experimental measurements are usually performed for a bundle of carbon nanotubes, if the 
tube-tube interaction is not weak it is also interesting to investigate the effect of inclusion 
of tube-tube interactions [29] . 
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FIGURE CAPTIONS 



Fig.l. The carbon-carbon two-body empirical potential as a function of the bond length 
for four sets of parameters. Our numerical analysis is based on the modified Brenner II 
potential. 

Fig. 2. Thermal conductance k versus temperature for different C-C bond lengths for zigzag 
(11,0) SWCN assuming perfect phonon transmission. Inset: The temperature-dependence 
of the thermal conductance for armchair (11,11) and zigzag (11,0). 

Fig. 3. Thermal conductance k versus temperature for different resonance width (D) fixing 
D = E Q in unit of cmT x . 

Fig. 4. Comparison of theoretical and experimental results. 
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